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Even before the experimental discovery of spin- and charge-stripe order in La2-c C -yNd y Sr cc Cu04 
and La2-ccBaccCu04 at x = 1/8, stripe formation was predicted from theoretical considerations. 
Nevertheless, a consistent description of the complex coexistence of stripe order with superconductivity 
has remained a challenge. Here we introduce a Hartree-Fock decoupling scheme which unifies previous 
approaches and allows for a detailed analysis of the competition between antiferromagnetism and 
superconductivity in real and momentum space. We identify two distinct parameter regimes, where 
spin-stripe order coexists with either one- or two-dimensional superconductivity; experiments on 
different striped cuprates are either compatible with the former or the latter regime. We argue that 
the cuprates at x = 1/8 fall into an intermediate coupling regime with a crossover to long-range 
phase coherence between individual superconducting stripes. 



PACS numbers: 74.72.-h,74.20.Rp,74.25.Ha 

I. INTRODUCTION 

Since the discovery of high-T c superconductivity in 
La2- x Ba :c Cu04 [1], the analysis of many experiments has 
led to a complex phase diagram of cuprate materials, in 
part originating from the competition between antifer- 
romagnetic (AF) and superconducting (SC) correlations. 
Both are induced by the Coulomb repulsion on the copper 
d-orbitals in the Cu02 planes. At a critical hole-doping 
level antiferromagnetism and superconductivity are in 
balance, and details of the material or its environment 
determine which order is realized or whether a regime 
of coexistence or local phase separation exists. At very 
low temperatures most cuprates show a transition from 
an AF to a SC state, if the density of charge carriers 
exceeds a critical value, and superconductivity vanishes 
again beyond a higher carrier density, when the Coulomb 
interaction becomes less significant. 

Above the critical hole doping level, AF correlations 
on reduced time and length scales persist and fluctuat- 
ing antiferromagnetism may coexist with Cooper pair- 
ing. One prominent phase with a combination of both 
has "stripe order" . An ordered phase with static stripes 
was first observed in 1994 by Tranquada et al. in the 
nickelate La2Ni04.i25 [2] and later on in the rare-earth 
doped cuprate La2-x-^Nd 2/ Sr :c Cu04 [3, 4] and in the orig- 
inal high-T c cuprate La2- x Ba x Cu04 [5-7]. The cuprate 
systems evolve from a Mott insulator into a homoge- 
neous SC phase as x is increased, but as x approaches 
1/8, antiferromagnetism returns in a characteristic spin- 
density- wave (SDW) pattern and superconductivity is 
suppressed [4, 8, 9]. Based on his neutron scattering data, 
Tranquada suggested the existence of antiferromagnet- 
ically ordered spin ladders separated by metallic lines, 
thus forming stripes with a width of four lattice constants 
(figure 1). The metallic lines form anti-phase domain- wall 
boundaries between the AF stripes, which optimizes the ki- 



netic energy through virtual electronic hopping processes 
into the AF stripes [10]. The appearance of static stripe 
order seems to be tied to the structural transition into 
the anisotropic low-temperature tetragonal phase (LTT) 
in which at low temperatures superconductivity possibly 
coexists with stripe order [9, 11, 12]. The spin dynamics 
of such systems, has been explained quite successfully us- 
ing models of coupled spin ladders [13-15], although their 
microscopic origin has remained unresolved within this 
ansatz, and superconductivity was not incorporated. In 
this framework both, site and bond centered (i.e. three- or 
two-legged) ladders are compatible with the experimental 
results for the dynamical spin susceptibility, although dis- 
tinct theoretical considerations suggest that site centered 
stripes are favored [15]. 

Unidirectional charge- and spin-stripe order was indeed 
predicted for cuprate systems before their experimen- 
tal discovery. In 1989, Zaanen and Gunnarson [16] and 
Machida [17] found AF stripe formation in the mean- field 
solution of a two-band Hubbard model. After stripes had 
been discovered experimentally, the t-J model was identi- 
fied as a more appropriate approach towards stripe forma- 
tion in hole doped cuprates. Extensive numerical studies 
were performed using exact diagonalization [18], DMRG 
methods [19, 20] or Monte Carlo simulations [21, 22], all 
indicating that stripes indeed form in the hole doped t-J 
model. Although these numerically exact results yielded 
a consistent picture, their range of validity is limited 
by small system sizes and boundary effects, and the in- 
clusion of superconductivity has only recently become 
possible [20, 22]. 

Meanwhile, new approaches have been developed for a 
phenomenological characterization of the SC state that 
possibly coexists with spin- and charge-stripe order. Berg 
et al. suggested a new type of SC state, termed "pair 
density wave" (PDW) , in which Cooper pairs with center- 
of-mass momenta q and — q coexist and the pair den- 
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Figure 1: Three-legged spin-ladder structure with an anti- 
phase domain wall as suggested by Tranquada etal. [3] for 
striped cuprates at hole doping x — 1/8. The ladders are half 
filled and antiferromagnetically ordered, whereas the domain 
walls are quarter filled and metallic. 



sity oscillates with wave vector q [23-27]. The latter 
oscillation is analogous to the FFLO solution of a su- 
perconductor in magnetic fields presented by Fulde and 
Ferrell [28] and by Larkin and Ovchinnikov [29], and is 
typically accompanied by a CDW with wave vector 2q. 
Houzet and Buzdin introduced a description of the FFLO 
state on the Ginzburg-Landau level [30] , which was later 
used as a model for the PDW by Agterberg and Tsunet- 
sugu [23] and by Berg et al. [25]. Subsequently, the PDW 
state without antiferromagnetism was shown to be the 
energetically stable solution for a microscopic nearest- 
neighbor pairing model with sufficiently strong pairing 
interaction [27]. This state indeed explains qualitatively 
some of the observed properties of striped cuprates [26]. 

To examine the nature of the SC state in coexis- 
tence with spin- and charge-stripe order, a simple mi- 
croscopic model is desirable which allows to adequately 
describe magnetic and SC order simultaneously, merging 
the strong-coupling ansatz of the t-J model with the weak- 
coupling framework of BCS theory. Work in this direction 
was recently performed using a Gutzwiller renormalized 
mean-field ansatz for the t-J model [31, 32]. However, 
the PDW did not prove to be the groundstate, but rather 
a "modulated d-wave" state without phase shift of the SC 
order parameter between neighboring stripes was found 
as the lowest energy solution. In contrast to the PDW, 
this "modulated <i-wave" has a pair density that oscillates 
around a finite uniform q = component, and pair and 
charge densities oscillate with the same wave vector. 

In this paper we introduce a modified t-J model as an 
extension of the standard t-J model which allows for a 
straightforward mean-field decoupling. This method of 
relaxing the no-double-occupancy constraint and extend- 
ing the t-J model by an additional on-site repulsion was 
originally introduced by Kagan and Rice [33] to demon- 
strate the possibility of d-wave superconductivity in the 
t-J model. This ansatz enables us to study stripe forma- 
tion and superconductivity in a wide range of parameters, 
and it merges two previously distinct approaches to the- 
oretically describe cuprates in different regimes of the 



renormalized exchange coupling J. 

In section II we introduce the t-J model and identify 
two parameter regimes where a mean-field description 
of magnetic correlations is appropriate for hole-doped 
systems: i) J <C t. In this limit superconductivity domi- 
nates and antiferromagnetism is controlled by the on-site 
Coulomb repulsion rather than by the superexchange in- 
teraction. For small J the t-J model therefore maps onto 
a simple nearest-neighbor pairing Hamiltonian with an 
additional on-site repulsion U ; this model is derived in 
section II B as the "[/-model", ii) J ~ t. A strong ex- 
change interaction tends to phase separate the system into 
a hole- free antiferromagnet and a hole-rich part [35-37]. 
We show that for J ~ t an intermediate regime exists in 
which the charge is separated periodically in half-filled AF 
regions and quarter- filled metallic regions. In this regime 
the Coulomb repulsion is less significant and the t-J 
model approaches the so called "V- model" , which was in- 
vestigated in detail in [34] in particular with respect to the 
breakdown of superconductivity in I^-^-yNdySr^CuC^ 
near x = 1/8. A similar regime has also been found in 
the t-J model itself using DMRG calculations [20] and 
variational Monte Carlo simulations [21]. 

Section III describes the numerical strategy to obtain 
striped solutions within the framework of the Bogoliubov 
- de Gennes (BdG) equations. Although charge- and spin- 
stripe order coexisting with superconductivity is found 
for J <C t as well as for J ~ t, there are important 
qualitative differences in the two regimes. We provide 
a detailed comparison for the origin of stripe formation 
in both regimes in section IV. The striped groundstate 
solution presented in this section is energetically favored 
among a considerable number of other possible solutions. 
Many of these energetically unfavorable solutions exist 
only in certain temperature ranges. We therefore present 
a survey of the temperature dependences of the most 
regular solutions in section V. In section VI we draw 
conclusions and discuss our findings in relation to striped 
cuprate materials. 



II. t-J VERSUS BCS MODEL 

Shortly after the discovery of high-T c -cuprates it was 
agreed upon that the electronic structure of the Cu02 
planes is well described by a three-band Hubbard (Emery) 
model [38] , which can be further reduced to a one-band 
repulsive Hubbard model for electron densities above 
2/3 [39]. Apart from numerically exact solutions, two 
different approaches have commonly been followed to 
treat the Hubbard interaction U^2 i s n^ s ni 5 _ s : i) In the 
strong-coupling regime (U ~ 10 1), the Hubbard model 
is approximated by the t-J model with a finite superex- 
change coupling J [40, 41]. This effective interaction is 
non-local and responsible for superconductivity as well as 
for antiferromagnetism. ii) Alternatively, effective non- 
local interactions in the finite U Hubbard model can be 
computed numerically using, for example, QMC Simula- 
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tions [42] or determined diagrammatically, e.g. in the 
RPA approximation. One thereby obtains a non-local at- 
tractive interaction between electrons on nearest-neighbor 
sites which may serve as a basis for a BCS type mean-field 
description. 

In this section, we derive mean-field Hamiltonians for 
both of the above approaches, i) and ii), and show that 
they indeed describe the same physics and are formally 
equivalent. 



A. t-J Model 



We start from the t-J Hamiltonian: 



n,j = - E tij ctc js + j E (s« • s, - ^) , (i) 



where n is = c\ s c is , m = n^+n^, and c M = <H i3 (l-ni- s ) 
is the projected annihilation operator which enforces the 
constraint that no lattice site is doubly occupied. The spin 
operator is S* = J2 s ,s> c \s°ss'Cis' /2 with u = (cr x , a y , cr z ) 
containing the Pauli matrices. Here we use the matrix el- 
ements tij = {t,t'} for nearest and next-nearest neighbor 
hopping with t' = — 0.4 £, and (i, j) denotes all pairs of 
nearest-neighbor sites i and j. Note that the three-site 
terms are omitted in the Hamiltonian (1). 

The constraint on the site occupation is the main obsta- 
cle in solving the t-J model and also prevents a straight- 
forward mean-field approximation. An ansatz to solve 
this problem is the slave-boson technique (see [43] for a 
review). Here we follow alternatively a method by Kagan 
and Rice [33] who introduced an extension of T-Lt-J, in 
which an additional Hubbard repulsion U replaces the 
local constraints: 

n t .j = - E ^ c L^> + jJ2( s i- s i - ^) 



the decoupled Hamiltonian becomes 

U 



t-j 



tij C l<s C 3,s + ^ ^1 n i,sCj^- s Ci,-s 
i,j,s i,s 

+ E[^t+^i ( 4 ) 



(i,3) 



^ (jrii- s c\ Q c 



with 



<^ijs = J{ c \ s C js)- 

The constant Ci is given by 



(5) 

(6) 
(7) 



Cv 



=E 

(id) L 



A S *A S - 

3113 



2 J + g ^t n j4 



/7 



E 



(8) 

Ci is essential for determining the free energy correctly. 
The superconducting order parameter A^- represents 
nearest-neighbor electron pairing in the spin-singlet chan- 
nel, while possible spin triplet contributions cancel within 
the Hartree-Fock decoupling of the t-J model. n^ s is the 
spin resolved thermal average of the local charge density 
rii and Xij S renormalizes the nearest-neighbor hopping 
amplitude t. 

The pairing term in W^j drives superconductivity with 
a maximum energy gain for an order parameter A^- with 
el-wave symmetry, i.e. A^±# = — A^±£. Furthermore, 
T~L^j contains two terms which favor antiferromagnetic 

order. The Hubbard term Uni jS c\ _ s q ? _ s costs the energy 

U for each doubly occupied lattice site and thus polarizes 
the spin on each site. The kinetic energy in turn is 
optimized by ordering the spins antiferromagnetically. 



U 



The term — Jrii 



C ],s C 3,s 



directly supports an opposite 



H y ^ nj^rii-s. (2) spin alignment on nearest-neighbor sites by the energy 



In the limit U — >• 00, one recovers Ht-J from (2) with 
J = J — 4t 2 /U, if the three-sites terms are omitted. Gen- 
erally, mean- field solutions of W, t _j overestimate antifer- 
romagnetism as compared to superconductivity at large 
interaction strength. Therefore more physical results are 
expected for smaller values of U (c.f. the discussion in 
section IV). 

The Hamiltonian T~L t _j is readily decoupled on the 
Hartree-Fock level. With 



Si- 



TliUj 



= - [ - n it n jl - n 4 n it + c^c}^ + cj^cj^ J 



(3) 



J\mi\(nj^ s + nj 5 _ s )/2 per bond, where ra$ = — 
is the local magnetic moment. For large values of J the 
system approaches perfect AF order and the Hubbard 
term looses relevance. On the other hand, even if all 
lattice sites are fully spin polarized by a large Z7, a finite 
J term in (4) leads to AF correlations. 

The above reasoning justifies to use l~L t _ j with a finite U 
instead of the original t-J Hamiltonian for the following 
physical reason: In the limit U — >• 00, the half-filled 
one-band Hubbard model decouples into independent 
localized electrons without any magnetic order. In the 
t-J model, this is reflected by J = t 2 /U —> 0. If a multi- 
band Hubbard model is mapped onto the t-J model, the 
superexchange terms lead to a finite and sizable J and 
thus to antiferromagnetic order even in the limit of a large 
intra-orbital U [40]. Antiferromagnetism is in this case 
controlled by J rather than by U. 
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B. U- and V-Mode\ 

The dominant effective interactions derived from the 
QMC calculations in [42] are the on-site repulsion U 
and an attractive interaction V for electrons on nearest- 
neighbor sites. If only these two dominant interactions are 
kept, we arrive at the real-space BCS type Hamiltonian 

t ^ t t 

*LuV — ~ / J ^ij c i^ s C hs ~ ~W / J C i,s C j,-s C j,-s c i,s 



U 

+ ~2 Z^ n i,sni -a, (9) 

i,s 

with V > and U > 0. The mean-field decoupling of 
T-Luv generates the bilinear Hamiltonian 



U 



tij C l<s C 3,s o n i,s c j,-s c i, 



UV 



with 

and 

C 2 = 



-In 



(ij) 



A* - A- • 



V 



+C 2 
(10) 

(11) 



u 



2 4^ 



(12) 



Comparing T~L t _j to Tiuv reveals that the mean-field 
decoupled Hamiltonians are almost identical, if one identi- 
fies the interaction parameters {J, U} < — > {V,U} in Ti^j 
and Hyy, respectively. In W§ v , the term X^ s is missing, 
which renormalizes the nearest-neighbor hopping t in the 
t-J model. Since it effectively decreases £, it tends to 
stabilize superconductivity; its qualitative influence on 
the groundstate solution is however marginal (c.f. [44]). 
More relevant are the differences in the definitions of 
the SC order parameters in Ti^j and T~L^ V : whereas A 3 - 
contains only the spin-singlet channel, includes also 
the S z = triplet channel. In a non-magnetic system, 
the triplet component generally vanishes for a nearest- 
neighbor interaction, because the pair wave-function has 
even parity. In the presence of antiferromagnetism how- 
ever, a finite triplet admixture appears (c.f. [44, 45]). 
In Hu V , this "7r-triplet" has its own order parameter 
A^- = -V((c^Q t ) + (c jt c^))/2, the spin-triplet compo- 
nent of Aij. However, A^- remains small compared to the 
spin-singlet component, as we have verified. Therefore we 
neglect the triplet component in the following calculations 
and discussions and use the spin-singlet order parameter 
A* ■ = —V((cjiCrf) — (ctfCii))/2 alone with the coupling 
constant V [instead of J as in (5)]. Furthermore, W^v has 



an additional prefactor 1/2 in the SC term. Therefore SC 
order is weighted stronger in W^j than in T€§ v and con- 
sequently J has to be larger than V in order to stabilize 
antiferromagnetism (see the discussion in Sec. IV A). 



Notice that the last term V ^Z s ^i,-s c 



'j, 8^3,8 



in (10), 



which appears identically also in Ti^j, is absent in the 
classical BCS theory, because in a homogeneous system it 
only renormalizes the chemical potential \i used to control 
the particle number (see section III). It is however a 
source of AF order in the nearest-neighbor pairing model 
T-Luv- 

A mean-field approach does not necessarily capture all 
the physics contained in the underlying model Hamiltoni- 
ans. The above decoupling scheme lacks e.g. the proper 
physics of the Mott insulator to metal transition. There 
are however two specific cases in which the mean-field 
solutions of T-Luv or ^t-J are m deed self-consistent: 

1. If V <C £, the term —Vfii- s c^ s Cj^ s contributes little to 
the emergence of antiferromagnetism, which in this case 
is controlled mainly by U. If also U is small, supercon- 
ductivity is the dominating order and the SC gap is larger 
than the AF gap. Therefore the failure of the mean-field 
theory to describe the quasiparticle peak at the Fermi 
energy of a doped Mott insulator is not relevant. Thus, 
for small V, W§ v reduces to the "[/-model" : 



u 



J2 n i,8c{_ 



+ o £ \ A UA + A ^ c ^tl (is) 



3,s^3,s 



discussed in detail in section IV A. 
2. If V is close to or larger than t, the term Vfii- 
is the dominant source of antiferromagnetism. The sys- 
tem separates into nearly half-filled AF regions (stripes) 
where superconductivity is suppressed, and into empty 
or quarter-filled metallic (or superconducting) regions, 
depending on the type of solution. Since commensurate 
antiferromagnetism appears only in the half-filled regions, 
a mean-field approach is consistent also in this case, and 
Tifj V reduces to the "V-model": 



Tiy — — y^tij c ts C 3\s 

i,3,s 



(i,3) 



(14) 



discussed in [34] and in section IV B. Within the V-model, 
a finite U -term can be added without significant influence 
on the groundstate solution, since the V-term alone is 
sufficient to drive the AF order parameter close to its 
maximum. 
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III. STRIPED SOLUTIONS 
A. The Bogoliubov - de Gennes Equations 

Here we briefly summarize the basic steps to solve the 
Bogoliubov - de Gennes (BdG) equations for and 
T~Ly and explain how striped (or other inhomogeneous) 
solutions are identified. 

The Bogoliubov transformation which diagonalizes the 
Hamiltonians (13) and (14), respectively, is given by 



n 



(15) 
(16) 



where the coefficients u n i S and v n i S are obtained from the 
eigenvalue equation 



A 

-it 



= En 



with u n , s = (u 

nis? • • • i u nis-> • • • ) &n.c! the corresponding 
vector v n , s . The sum over n in (15) and (16) extends 
over all positive eigenvalues E n and the corresponding 
eigenvectors. The operators t s and A act on u n s and 

v n,s as 



where j labels all nearest-neighbor sites of i and / labels 
all sites for which tu is finite, i.e. nearest and next-nearest 
neighbors. 

Using the symmetry of the energy spectrum E n around 
E = 0, it is sufficient to diagonalize (17) for a single spin 
component, say s = t- Inserting the transformation (16) 
into (6) and (11) leads to the self-consistency equations 
for the order parameter and the spin resolved densities 
n it and 



n 

(20) 

% = XI U niff( E n ~ M), (21) 
n 

n^l-J^/^-/*). (22) 

n 

Here we used that J2 n v^f(E n — fi) equals the number 
of holes with spin s = |, and f(E — fi) denotes the 
Fermi function. The chemical potential \i fixes the total 
particle number N to the desired value. Note that in 



the presence of antiferromagnetism, A* ■ has a finite spin- 
triplet component. 

The BdG equations (17) and (20 - 22) are solved iter- 
atively until self-consistency in A* • and fii s is achieved. 
In each loop fi is adjusted to keep TV constant. To ob- 
tain a higher momentum and energy resolution, we use 
the super cell method to block-diagonalize the eigenvalue 
problem (17). This procedure was introduced by Wang 
and MacDonald [46] for magnetic super cells (for more 
details, see [44] and [47]). Here the size of a super cell 
must be chosen commensurate with the wavelength of the 
striped solution we attempt to find. 

The BdG equations typically have more than one solu- 
tion, into which the self-consistency cycle may converge. 
An anticipated solution can usually be selected by choos- 
ing appropriate initial values for A* • and fii S (and for /i), 
the search for new solutions is however always a demand- 
ing task. Typically, two types of SC solutions compete in 
the presence of AF stripes: the "modulated d-wave" and 
the PDW solution. Both are characterized by an order 
(17) parameter with local d-wave symmetry, i.e. A* i±& has 
the opposite sign of A* but the absolute values are 
different, if an extended s- wave component exists. In the 
PDW solution, A^- also has opposite signs on neighboring 
stripes. Such sign changes in A^ or in the antiferro- 
magnetic order parameter Mi = (— l) z ra$ do not evolve 
continuously from a uniform initial state. In order to trace 
the anti phase striped solutions discussed in section IV, 
the initial values of A^ and Mi must therefore have the 
same anti phase stripe pattern with equal wavelength. 

The groundstate for each set of parameters is found by 
minimizing the system's free energy F = (H) — TS over 
all self-consistent solutions. This means minimizing F 
over stripe patterns with different wavelength, but also 
over bond- and site-centered stripes and different types 
of SC states, i.e. a "modulated d-wave" or a PDW state. 
Three general observations are made: 

1. For those values of V and U for which striped 
solutions exist at all, solutions of the BdG equations exist 
for all stripe widths that are commensurate with the 
finite-size lattice, with AF stripes separated by an anti 
phase hole-rich domain wall. 

2. For x close to 1/8 the striped groundstate has a 
stripe wavelength A = 8a, where a is the lattice constant. 
For small values of V and U, the groundstate is a homo- 
geneous non-magnetic d-wave superconductor, whereas 
phase separation occurs, if V or U exceeds a certain limit. 

3. A SC state of PDW type occurs only, if the pairing 
interaction V is large enough (V > t). Its energy however 
is always larger than the energy of a "modulated d-wave" 
state, although the energy difference becomes vanishingly 
small for large V. This result is in agreement to ear- 
lier work using either a Gutzwiller approximation [32] or 
DMRG calculations [20] for the t-J model. In contrast, 
the pure PDW in non-magnetic systems was shown to 
be the SC groundstate for a sufficiently strong nearest- 
neighbor pairing interaction [27]. 

A thorough analysis of the groundstate properties in 
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Figure 2: Real-space characterization of the "modulated d-wave" (mdSC), PDW, and non-SC solutions of the [/-model, (a) Local 
magnetization rrn = — n^. (b) Charge density m — + n^. (c) d-wave projection A^ of the SC order parameter, (a - c) 
correspond to the PDW solution, (d) AF order parameter Mi — ( — lYrrii. (e) Charge density rii. (f) SC bond order-parameter 
All results were obtained on a 16x12 lattice for V — 2t, and U — 3.4 1. 



different parameter regimes follows in section IV. All 
calculations, except for section V, were performed at 
temperature k^T = 0.01 1. 

B. Comment on Energy Minimization 



a variable particle number N. However, our calculations 
aim at solutions with fixed N by adjusting /i, although all 
expectation values are evaluated grand canonically. There- 
fore the thermodynamically stable state is determined by 
the minimum of F = + fiN . 



In the grand canonical framework of BCS-type mean- 
field theories, the thermodynamically stable groundstate 
is defined by the global minimum of the grand canonical 
potential f2 = (H) — TS — jiN where \i is the chemical 
potential and the entropy S of the system is given by 

S = -k B ^2 if ^n) In f(E n ) + f(-E n ) In f(-E n )} . 

n 

(23) 

Self-consistent solutions of the BdG equations correspond 
to local minima in Q. To determine the groundstate or the 
thermodynamically stable state at finite temperatures, it 
is therefore necessary to compare different self-consistent 
solutions of the BdG equations for each chosen set of 
parameters (e.g., with or without sign changes in A^-, 
site-centered or bond-centered stripe patterns, etc.), and 
to select the solution with the lowest Q. 

This minimization procedure is valid for the grand- 
canonical ensemble with a fixed chemical potential \i and 



It is instructive to compare Q and F for solutions of the 
BdG equations with fixed \i or with fixed N, respectively. 
For fixed N, one can prepare the system initially in a state 
where Q is minimal, but not F. Specifically this would 
be a state where antiferromagnetism is absent. For fixed 
TV, Q has a global minimum without AF order, whereas 
F has two degenerate minima with AF order. As F is 
not minimal for this state, the non-magnetic solution is 
unstable and the BdG equations will eventually converge 
into one of the two minima of F. By contrast, the roles 
of ft and F are interchanged, if ji is fixed and N varied. 
Therefore fixing TV instead of \i is numerically equivalent 
to a Legendre transformation from the grand canonical 
back to the canonical ensemble. 



Figure 3: Momentum-space characterization of the "modulated 
d-wave" solution of the U model for x = 1/8, V = 2 1, and U = 
3 At on a 16x12 lattice with 7x7 supercells. (a) Occupation 
probability n(k), (b) pair density -P(k), (c) spin order ps(k), 
and (d) charge order pc(k). 



IV. STRIPE PATTERNS IN THE U- AND THE 
y-MODEL 

A. [/-Model 

If V <C t (say V ~ 0.2 1), the energy gain originating 
from — V ^2 S n^-gC^^Cj^ upon ordering antiferromagnet- 
ically is smaller than the accompanying cost in kinetic 
energy. If also U < t, the groundstate will be a homo- 
geneous, non-magnetic d-wave superconductor, whereas 
larger values of U induce AF order. For moderate values 
of U, d-wave superconductivity is still the dominating or- 
der, which coexists with weak antiferromagnetism above 
a critical value of U depending on the pairing interaction 
strength. This regime is described by the [/-model given 
by the simplified mean-field Hamiltonian in (13). 

The L/-model supplemented with non. magnetic impu- 
rity potentials has been widely and successfully used 
to describe AF correlations in disordered cuprate super- 
conductors [47, 49] or around vortex cores in magnetic 
fields [44, 46-48]; Andersen and Hedegard showed the 
existence of striped solutions for this model [50] . In partic- 
ular disorder induced antiferromagnetism appears above 
a critical value of U. Here we show that this model also 
allows for striped groundstate solutions in clean systems 
for U c i < U < Z7 C 2, where the critical value U c \ ~ 3V 
is slightly larger than for anitferromagnetism in disor- 



Figure 4: The same quantities as in figure 3 but for the PDW 
solution for the same parameter set. 



dered systems. For values of U > U C 2 ~ 6 V, we obtain 
a homogeneous AF solution without superconductivity, 
i.e. a doped AF Mott insulator, for which the mean-field 
treatment is not adequate. 

Figure 2 shows the typical groundstate solution for 
hole doping x = 1/8, V = 2t, and U = 3 At. The 
local magnetization = — (a) forms the well 
known site-centered spin-stripe structure with an anti- 
phase domain wall between the stripes and an overall 
periodicity of 8 a. Upon increasing U through f7 c i, the 
maximum polarization max^ \rrii\ abruptly increases to 
^0.5, and rises further to ~0.75 towards C/ C 2- Thus the 
system is quite far from the ideal spin-ladder structure 
with fully polarized AF stripes (in contrast to the V- 
model, c.f. section IV B). Holes are expelled from the AF 
stripes, leading to a charge density modulation with a 
relative amplitude of 5-10% [see figure 2 (b) and (e)]. 

A solution with bond-centered stripes in a two-legged 
spin-ladder structure also exists, but it is higher in en- 
ergy for the parameters chosen above. If however U is 
increased towards f7 C 2? the groundstate changes to the 
bond-centered spin structure. 

The groundstate for U c \ < U < U C 2 is a "modulated 
d-wave" (mdSC) state. The d-wave projection of the SC 
order parameter 

[see figure 2 (c)] shows a similar stripe pattern as the 
charge density. Because of the x - y asymmetry of the 
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Figure 5: Momentum-space characterization of the non-SC so- 
lution of the U model for x = 1/8, V = 2t, and U — 3.4 1 on a 
16x12 lattice with 7x7 supercells. (a) Occupation probability 
n(k), (b) spectral function A(k,ou = 0). 



striped system, a finite extended s-wave component 

A* = (A% i+& + A% i+i) + A' M _ 4 + (25) 

is also induced. To fully characterize the SC state, x and 
y bond-order parameters are shown in figure 2 (f) over a 
cross section of two stripe periods. For comparison, we 
also include a solution of the PDW type with a sign change 
in Af • from one AF stripe to the next. The PDW solution 
has a reduced SC order parameter as compared to the 
mdSC state and a slightly higher free energy. Although 
the maximum AF order parameter Mi is larger for the 
PDW solution, the gain in magnetic energy is not sufficient 
to compensate for the smaller SC condensation energy 
caused by the zeros in the SC order parameter in the 
center of the AF stripes [see figure 2 (f)]. 

If superconductivity is completely suppressed (by set- 
ting Afj = artificially), the modulations in and 
rti are barely larger than in the mdSC state and almost 
identical to the PDW solution [see figure 2 (d) and (e)]. 
This indicates that the spin polarization is limited mostly 
by the kinetic energy, which energetically favors a ho- 
mogeneous charge and spin density, and less limited by 
the competition with superconductivity. Indeed, anti- 
ferromagnetism and superconductivity rather avoid each 
other, as becomes evident from the momentum-space 
characterization discussed below. 

A complementary description of the intertwined SC 
and AF orders, together with the resulting anisotropy 
of the system, is possible in momentum space. Note 
that although we measure all length in real space in 
units of the lattice constant a, we will set a = 1 for 
the following discussions in momentum space. Figure 3 
characterizes the groundstate solution of the mdSC state. 
The occupation probability n(k) in figure 3 (a) and the 
pair density P(k) in figure 3 (b) are almost identical to 
the known distributions of a pure d-wscve state without 
magnetism (c.f. [27]), except for a small x-y asymmetry. 



The pair density squared 

P2 ( k ) = [( C kt C kt C -k+q| C -k+q|) 

q 

- (4t c kt) (c f _ k+qi c_ k+qi )] (26) 

measures the average correlated occupation of the elec- 
tron states with momenta k and — k + q. Thus P 2 (k) 
is largest at those momenta in the Brillouin zone where 
electron pairs predominately form. Within the mean-field 
decoupling scheme, equation (26) reduces to P 2 (k) = 
Xlq( c -k+q|Ckt) 2 - The pair density is concentrated 
around the Fermi surface of the normal state and vanishes 
at four nodal points near the zone diagonal. 

Around these nodal points, AF correlations emerge. 
They are quantified by the spin order in momentum 
space, that is by p s (k) = E q E s ( 5C k+ qs c ks) [ see fig- 
ure 3 (c)]. ps indeed has its maxima at the nodal points, 
reflecting the competition between, or rather avoidance 
of superconductivity and antiferromagnetism, i.e. anti- 
ferromagnetism is strong where pair formation is weak. 
Pc(k) = Eq/o E s ( c k+ qs c k s ) [see figure 3 (d)] is the ana- 
log for the charge order, and pc(k) is non zero only in 
the presence of charge modulations. 

It is instructive to compare figure 3 to the corre- 
sponding results for the PDW solution in figure 4 or 
to a state where superconductivity is suppressed arti- 
ficially (figure 5). If superconductivity is absent, the 
spin polarization is slightly stronger in the AF stripes 
which leads to the characteristic occupation probabil- 
ity n(k), shown in figure 5 (a): most of the Fermi sur- 
face is gapped by the AF order with a continuous n(k) 
across, but discontinuities in n(k) remain on the border 
of a horizontal bar between k y = 7r/4 and k y = — 7r/4. 
These discontinuities constitute disconnected Fermi sur- 
face arcs which are clearly visible in the spectral function 
A(k,uj = 0) = -(l/7r)ImG(k,k,cj = 0) shown in fig- 
ure 5 (b). 

In the superconducting PDW solution (see figure 4), 
antiferromagnetism is almost as strong as in the absence 
of superconductivity [c.f. figure 2 (d)], thus n(k), shown 
in figure 4 (a), has the same structure as in figure 5 (a), 
but the entire Fermi surface is now gapped due to super- 
conductivity. The pair density P(k) of the PDW state 
shown in figure 4 (b) is strongly anisotropic because of 
the line zeros of the order parameter enforced by the tt 
phase shift between neighboring stripes. -P(k) is similar 
to the pure PDW without antiferromagntism described 
in [27], the latter solution is however not gapless. The 
pairing free pieces of the reconstructed Fermi surface of 
the pure PDW [24, 26, 27] are now gapped due to anti- 
ferromagnetism and extend towards the (0, ±7r) points 
where superconductivity is weak for vertically oriented 
stripes [see figure 4 (c)]. Although the pairing contribu- 
tions around (0, ±7r) are weak compared to the vicinity 
of (±7r, 0), they are essential for a globally phase coherent 
SC state. 
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Figure 6: Pairing and effective hopping amplitudes in real 
space for the U model and the same parameters as in figure 4. 
(a) Pairing amplitude P(r^) = (co^c^). The value of the on- 
site term is P(0) = 1.88. (b) Effective hopping amplitude 
pc(?i) — ^2 s ( c 0s c is)- The value of the on-site term is pc(fi) = 
1.73. 



It is instructive to analyze also charge and pair corre- 
lations in real space. The effective hopping amplitude 
pc( r i) = ^Zs( c ^s c is) f° r hopping from site to site i is 
obtained by Fourier transforming n(k) + pcQs). The pair- 
ing amplitude in real space is P(r^) = (cq^c^). Figure 6 
shows P{vi) and pc{ Y i) fc> r the /7-model. Both quanti- 
ties exhibit an oscillating behavior in which P(r$) decays 
over the distance of the SC coherence length which 
is roughly 10 lattice constants for the chosen parameters. 
The oscillations in Pfa), as well as the Friedel type of 
oscillations in pc(r^), extend over almost the same length 
in x and y direction, which verifies the 2D character of 
the solution of the [/-model. 

The dominance of the SC order in the /7-model is most 
apparent in the local density of states (LDOS) on mag- 
netic and non-magnetic lattice sites as shown in figure 7 
for the same parameters as above. For the "modulated 
d-wave" solution [see figure 7 (a)] the LDOS vanishes 
linearly upon approaching the Fermi energy, which is the 
characteristic of <i-wave superconductivity, but it has a 
slightly asymmetric shape, caused by the particle-hole 
asymmetry induced by the presence of a charge density 
modulation. The LDOS of the PDW solution [figure 7 (b)] 
has a smaller but s- wave like gap, since it has a strong 
extended s-wave contribution. In both solutions the gap 
scales with the SC order parameter, and its site indepen- 
dence is characteristic for the /7-model. If U is further 



Figure 7: LDOS for the U model and the same parameters as 
in figure 4. (a) mdSC solution, (b) PDW solution. 



increased to values above Z7 c2 , the magnetism induced gap 
starts to dominate over the SC gap and superconductivity 
breaks down. The /7-model therefore does only allow for 
SC solutions with superconductivity dominating over AF 
order. This behavior is typical for the /7-model and quite 
different to the y-model, as we show below. 



B. V-Mode\ 



A detailed description of the V-model Hy at x 



1/8 

was given in [34] . Here we extend this analysis to clarify its 
relation to the /7-model and the characteristic differences. 
In the F-model, the interaction strength V controls both 
the SC and the AF correlations. Upon turning on and 
increasing V, we find a similar sequence of phases as in 
the /7-model with increasing U at x = 1/8: For small V, 
the groundstate is a homogeneous, non-magnetic d-wave 
superconductor. At a critical value V c \ ~ 0.8 1, AF stripe 
order with wavelength 8 a sets in, similar to the solution 
of the U model; this solution is characterized in real space 
in figure 8forV = 2tatx = 1/8. Above a second critical 
value V C 2 ~ 2.6 1, superconductivity disappears and the 
system becomes susceptible to phase separation. 

Figure 8 shows that the AF stripes are almost half 
filled with a staggered spin polarization that is close 
to maximal, whereas superconductivity vanishes in the 
center of the AF stripes for both the x and y bonds. 
The system can therefore be described as locally phase 
separated into half- filled AF stripes and ID metallic lines 
acting as anti phase domain walls. The crossover from 



10 



AF order parameter Mi 



Charge density rii 



SC bond order parameter A* 




t y-b< 


)nds 






A 


A 


A 


A 


IX 


ft 


ft 


ft 




x- bonds 







(c) 

Figure 8: Real space characterization of the "modulated d-wave" solution of the V- model for x = 1/8 and V = 2 1 on a 16 x 12 
lattice, (a) AF order parameter Mi = (b) charge density m, and (c) SC bond order parameters A s ^ +ce and A s i;i+ ^. 




Figure 9: Momentum-space characterization of the non-SC 
solution of the V model for x — 1/8 and V = 2t on a 16x12 
lattice with 7x7 supercells. (a) Momentum occupation proba- 
bility n(k), (b) spectral function A(k,cj = 0). 



the homogeneous d-wave superconductor to the striped 
solution with amplitudes as in figure 8 is rather sharp. 

The ID character of the V-model solution is most ap- 
parent in momentum space, if superconductivity is again 
artificially suppressed. The momentum distribution n(k) 
[see figure 9 (a)] consists of a horizontal bar with occu- 
pied states for — 7r/4 < k y < 7r/4, and a diffuse cloud 
of occupied states in the center of the Brillouin zone far 
below Ep. The states with momenta inside the bar arise 
from the quasi ID metallic lines between the AF stripes, 
whereas the cloud near the zone center originates from 
the AF stripes themselves. Both structures are tied to 
separate energy windows in the LDOS [c.f. figure 11 (a)]. 
The borders of the bar form two ID Fermi surfaces at 
k y = ±7r/4, visible in the spectral function A(k, uo = 0) 
in figure 9 (b). Pairing can occur exclusively around 
these two Fermi surfaces and is therefore restricted to the 
metallic hole-rich lines. 

As for the [/-model, we also calculated the pairing 
and effective hopping amplitudes P(r^) and pc( r i) m 
real space, shown in figure 10. The pairing amplitude 
P(Ti) exhibits a SC correlation length which is highly 
anisotropic with £q « 2 a. The Cooper pair motion is 
therefore restricted to ID channels along the stripes and 
consequently the phases of the SC order parameters on 



adjacent non-magnetic lines become decoupled. In other 
words, the 2D phase coherent state is degenerate with un- 
correlated ID superconducting stripes. A globally phase 
coherent state may still be favorable by a Josephson cou- 
pling of the ID SC stripes, however on the Hartree-Fock 
level Josephson coupling is not included. Since charge 
transport in the x-direction transverse to the stripes is 
absent, also the charge correlations do not extend from 
one stripe to the next, as is visible from pc( r i) m fig- 
ure 10 (b). The dimensionality of the SC state therefore 
constitutes a qualitative disparity between the U- and 
V-model. 

The formation of energetically separated states on 
metallic (or SC) and AF sites is the essential distinc- 
tion between the V model and the U model. The LDOS 
of the non-SC system in figure 11 (a) shows a broad 
metallic band on the metallic sites around Ep, whereas 
on the AF sites all states are far from Ep. A SC gap 
therefore appears only on the non-magnetic sites [see fig- 
ure 11 (b)]. If V is lowered below V c i, the magnetism 
induced gap, which reaches from ~ — 2t to t in figure 11, 
becomes smaller than the SC gap, as in the LT-model. 
Here however magnetism vanishes below V c \. The V- 
model behaves therefore contrarily to the /7-model, where 
superconductivity is the dominant order until it vanishes 
when the AF gap becomes too large. 



V. TEMPERATURE DEPENDENCE 

In section III we mentioned already the existence of 
different self-consistent solutions of the BdG equations 
in both the U- and the V- model. They correspond to 
different local minima in the parameter space of the free 
energy F above the global minimum representing the 
states described above. Some of these additional solutions 
can only be obtained in special parameter ranges, or 
within a specific temperature regime. In this section we 
present the temperature evolution of a selection of these 
states. The list of states presented here is incomplete, but 
contains the most stable and regular patterns that are 
obtainable in both the U- and the V-model 

Whereas in the V-model the contribution to the free 
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Figure 10: Pairing and effective hopping amplitudes for the 
V model and the same parameters as in figure 9. (a) Pairing 
amplitude P(r^) = (co^c^), (b) effective hopping amplitude 
pc(ji) — ^2 s { c os c is)- The value of the on-site term is pc(0) = 
1.58. 



energy F from magnetism is by far larger than the conden- 
sation energy of superconductivity, in the [/-model both 
contributions are of similar size in the range of /7-values 
where superconductivity and AF coexist. Therefore local 
minima in F with distinct magnetic structures are much 
closer in energy and thus a large variety of solutions is ob- 
tained. On the other hand, in the V-model the dominant 
AF order allows only for a few self-consistent solutions. 

Figure 12 shows the maximum values of the SC or- 
der parameter A* • (a), the AF order parameter Mi (b), 
and the maximum amplitude of the CDW Srti (c) as a 
function of temperature for the /7-model for x = 1/8. 
For better characterization of the various solutions of 
the [/-model, we present the corresponding real-space 
patterns of ni and A^ in the upper panel. For the 
U -model we find solutions with stripe order of wavelength 
A = 8 a and A = 6 a, and also homogeneous solutions. The 
striped solutions exist both in the PDW and the mdSC 
type, whereas homogeneous antiferromagnetism may co- 
exist with homogeneous d-wave superconductivity (dSC). 
Some solutions disappear abruptly at certain tempera- 
tures. More specifically, superconductivity vanishes at 
temperatures where antiferromagnetism still persists. It is 
therefore instructive to trace the high-temperature states 
without superconductivity back to lower temperatures. 
These special solutions are indicated by white fields in the 
upper panel of figure 12 where the corresponding order 
(superconductivity or antiferromagnetism) is suppressed. 



Figure 11: LDOS for the V model and the same parameters 
as in figure 9. (a) Non-SC state and (b) SC state. 



Within the /7-model, the order for the onset of super- 
conductivity and antiferromagnetism is not fixed. Close 
to U c i, superconductivity sets in at a higher temperature 
than antiferromagnetism, while for the parameters used in 
figure 12, antiferromagnetism survives to slightly higher 
temperatures. In this regime, the sudden disappearance 
of superconductivity is characteristic for the [/-model, 
because, as discussed above, superconductivity breaks 
down, if the SC gap becomes smaller than the AF gap 
as the temperature is raised. This phenomenon is visible 
for the A = 8 a AF + PDW and AF + mdSC solutions 
as well as for the A = 6 a AF + mdSC solution. As the 
mdSC state vanishes with raising the temperature, the 
AF stripe order changes spontaneously from a site- to a 
bond-centered pattern. A remarkable observation is the 
disappearance of antiferromagnetism in the A = 6 a PDW 
solution. This state is stable only in a certain temperature 
window, while for lower temperatures it decays into an 
AF tartan pattern. 

For the V-model the magnetic order parameter is dom- 
inant and at x = 1/8 only solutions with wavelengths 
of 8 a or 16 a can be stabilized. The temperature de- 
pendence of these solutions is shown in figure 13. Here 
we also find non-superconducting states with AF stripes 
(A = 8 a SDW and A = 16 a SDW) as well as the striped 
superconducting state with wavelength A = 8a (A = 8a 
SDW+PDW); note that the phase difference of A^ on 
neighboring stripes is not uniquely determined in the 
mean-field solution of the V-model as discussed in detail 
in [34]. Since the non-magnetic sites in the A = 16 a SDW 
state are nearly empty, superconductivity is absent in 
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Figure 12: Upper panel: Selection of (meta) stable solutions of the [/-model. The groundstate is the A = 8 a AF + mdSC 
solution. In white fields, the corresponding order is excluded from this solution. Lower panel: Temperature dependence of 
the maximum SC order parameter (a), the maximum local magnetization (b), and the maximum amplitude of the CDW (c), 
characterizing stable solutions of the [/-model with U — 4 At and V — 2t. Some solutions are stable only in specific temperature 
windows and disappear where the corresponding lines end. 



this solution. A A = 6 a SDW state can be stabilized for 
higher doping levels, e.g. x = 1/6. Within the V- model, 
the superconducting T c is typically much smaller than the 
onset temperature of antiferromagnetism, and the two 
temperatures become equal, if V is lowered to V c \ . The 
A = 8 a SDW+PDW solution merges smoothly into the 
non-SC A = 8a SDW solution above T c , since Mi and 
Srii are almost temperature independent at T c . 



The sequence of phase transitions with decreasing tem- 
perature is complicated in the striped cuprate materials. 
It is best investigated in La 15 / 8 Ba 1 / 8 Cu04 [12, 52], where 
a CDW sets in at Tcdw ~ 54 K at the transition into the 
LTT phase, which is joined by an AF stripe order at a 
slightly lower temperature. Superconductivity appears 
only below 10 K, which suggests the dominance of AF 
order as in the F-model. A direct comparison to our 
mean- field models is however difficult, since these models 
do not describe fluctuating superconductivity, e.g. non 
phase-coherent stripes, which is likely to be present at 
much higher temperature close to Tcdw- 



VI. CONCLUSIONS 

In this article, we have explored how far a mean- field ap- 
proach of the t-J model and of closely related models can 
assist our understanding of static spin- and charge-stripe 
order observed in underdoped cuprates in coexistence 
with superconductivity. For hole doping x = 1/8 the free 
energy is minimized by an SDW order with wavelength 
A = 8 a with concomitant CDW order with A = 4 a, which 
agrees with neutron and X-ray scattering experiments for, 
e.g., La 15 / 8 Ba 1 / 8 Cu04 or La 15 / 8 _ 2/ Nd 2y Sr 1 / 8 Cu04. This 
type of order is stable for a wide range of the interaction 
parameter J in the t-J model and in both its descendants 
for small and large J, the /7-model and the V- model, 
respectively. These latter models represent two limiting 
cases where a mean-field treatment of the t-J model is 
well justified. 

While the energetically favored spin structure is a ro- 
bust property of the 7/8 filled t-J model, the nature of the 
SC state coexisting with the spin order changes between 
different limits. For small values of J, an almost homo- 
geneous and isotropic d-wave SC state is realized which 
recedes to one dimensional non-magnetic lines for larger 
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Figure 13: Temperature dependence of the maximum SC order parameter (a), the maximum local magnetization (b), and the 
maximum amplitude of the CDW (c), characterizing stable solutions of the V- model with V = 2t. 



J and also in the V-model. For larger J, phase coherence 
perpendicular to the stripes is lost or very fragile, which 
naturally explains the breakdown of superconductivity in 
the presence of static spin-stripe order, stabilized by the 
anisotropic hopping amplitudes in the LTT phase of rare 
earth doped 214 cuprates [4, 8, 32, 51]. If phase coherence 
is established with decreasing J, the favored SC order 
parameter does not acquire a phase shift between neigh- 
boring stripes and therefore has a finite homogeneous 
q = component. It is therefore not of the pure PDW 
type as suggested in [25], but rather similar to the state 
found in [32]. A definitive answer to the phase relation 
of the superconducting order parameter on neighboring 
stripes cannot be given within the presented mean-field 
description. A phase sensitive Josephson-coupling term 
is needed, in the spirit of the phenomenological theory 
of [25]. 

The analysis of the t—J model led us to two qualita- 
tively distinct regimes for the coupling constant J. Its 
appropriate value for high T c cuprates is generally be- 
lieved to be in the range of ~ 0.3 - 0.5 1. The physics 
of the cuprates therefore falls somewhere in between 
the two limiting regimes where the two mean-field ap- 
proaches are reliably applicable. While the two regimes 
are not continuously connected on the mean-field level, 
the model systems may well show a smooth crossover 
and features of both regimes. For example, ARPES mea- 
surements on La 15 / 8 Ba 1 / 8 Cu04 [53] show four-fold sym- 
metric Fermi-arc like structures at zero energy, compati- 
ble with the "almost" isotropic d-wave solution, whereas 
La 15 / 8 _ 2/ Nd 2y Sr 1 / 8 Cu04 with y = 0.6, which has a much 
stronger anisotropy in the hopping amplitudes along the 



x- and the ^-direction, shows a quasi one dimensional 
Fermi surface with a breach at the Brillouin-zone cen- 
ter [54]. These features fit well to the solution of the 
V-model in [34]. 

A related issue is the influence of disorder and vortices 
on the phenomenon of stripe formation. Our initial calcu- 
lations for disordered systems indicate that impurities act 
as pinning centers for AF stripes by slowing down spin 
fluctuations and thereby freezing static stripe order. In 
the £7-model, the stripe structures flexibly adjust in wavy 
forms to the impurity pattern, while superconductivity 
seems only slightly affected. On the other hand, impuri- 
ties with a moderate scattering potential do not disorder 
the stripe pattern in the V-model at all, although they 
are detrimental for superconductivity. Since stripes in the 
cuprates are not perfectly ordered, the V-model limit is 
probably not reached, but rather an intermediate regime. 
In the presence of antiferromagnetic stripes, the motion 
of vortices is confined to the stripe direction. Their stable 
position is at the center of an antiferromagnetic stripe, 
where superconductivity is weakest. As has been shown 
within the /7-model for homogeneous d-wave supercon- 
ductors [46-48], antiferromagnetism is enhanced around 
the vortex core, which, similarly to a strong impurity 
potential, forces the neighboring superconducting stripes 
to bend around the vortex. 
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